Entanglement of macroscopic test masses and 
the Standard Quantum Limit in laser interferometry 
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We show that the generation of entanglement of two heavily macroscopic mirrors is feasible with 
state of the art techniques of high-precision laser interferometry. The basis of such a demonstration 
would be a Michelson interferometer with suspended mirrors and simultaneous homodyne detections 
at both interferometer output ports. We present the connection between the generation of entangle- 
ment and the Standard Quantum Limit (SQL) for a free mass. The SQL is a well-known reference 
limit in operating interferometers for gravitational-wave detection and provides a measure of when 
macroscopic entanglement can be observed in the presence of realistic decoherence processes. 

PACS numbers: 42.50.Xa, 42.50.Lc, 03.65.Ta, 03.67.Mn 



The continuous quantum measurement of macro- 
scopic objects was first investigated in the context of 
gravitational- wave (GW) detection Q. In laser interfer- 
ometer GW detectors, incoming GWs induce very weak, 
yet highly classical, tidal forces on mirror-endowed test 
masses. These are approximately free masses since they 
are hung from seismic isolation stacks as pendulums 
with eigcnfrequcncy 1 Hz) much below detection band 
(lOHz < n/{2TT) < lOkHz). The interferometer mea- 
sures the change in the test-mass mirror's relative po- 
sitions 0]. The position observable, however, does not 
have commuting Heisenberg operators at different times, 
resulting in a Standard Quantum Limit (SQL) [l[ 
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where m is the reduced mirror mass. After a rather con- 
troversial debate it was eventually realized that the SQL 
can in principle be surpassed 0] if (i) measurement of the 
external force is made through monitoring of a quantum 
non- demolition observable (ii) the quantum mechan- 
ics of the measuring device is taken into account, and ap- 
propriate quantum correlations are used, e.g., with back- 
action evasion techniques To reach the sensitivity 
required for GW astronomy, future GW detectors as well 
as near term prototypes will be designed to reach and 
even surpass the SQL @, 0, H, 0, using either of the 
approaches mentioned. At the same time it is necessary 
to reduce all technical noise in order to reach this goal. 
These efforts are made with devices having test masses 
and arms ranging from CTam-scale and meter-long ho| , 
kilogram-scale and 100 m Q up to 40 kg and 4 km 

In this Letter, we show that a Michelson interferome- 
ter at its SQL produces entanglement between the con- 
ditional states of its otherwise free test-mass mirrors, if 
both common- and differential-mode test-mass motion 
are measured, but with different quantum-measurement 
processes - a direct mechanical analog y to the Einstein- 
Podolsky- Rosen gedanken experiment ll|. Contrary to 



former proposals [l3|, our work is aimed directly at ex- 
perimental efforts in the GW community towards SQL 
sensitivity, and therefore may enable the first observa- 
tion of entanglement between truly macroscopic objects. 
In technical terms, we introduce Wiener Filtering [isj . 
which, compared with the more conventiona l ap proach 
involving Stochastic Master Equation (SME) [1J|, allows 
us to determine the conditional state much more directly 
from the interferometer's noise budget. 

test-mass mirror (north) 



diff. mode 

photo detection 
com. mode (west) 



laser input 



Faraday 
rotator 



\3 



I eom. mode 



com. mode 



, test-mass 
' mirror (east) 



|iholo detection 
diff. mode (south) 



diff mode 



FIG. 1: Schematic plot of a power-recycled (PR) Michelson 
interferometer. Suspended test-mass mirrors are much lighter 
than other suspended optics. Differential motion of test-mass 
mirrors detected at dark (south) port and common motion at 
bright (west) port. A Faraday rotator might be used to access 
all of the back reflected light. 

Configuration. We consider an equal-arm Michelson in- 
terferometer as in Fig.m with laser (carrier) light injected 
from left. The two beams split at the beam splitter (BS) 
are reflected by identical mirrors in the north (n) and east 
(e) arms, before being recombined at the BS. The south 
port is kept dark at the zero point, with all light reflected 
to the bright port. A power-recycling (PR) mirror is po- 
sitioned such that it forms a resonant, but relatively low 
finesse, cavity for the carrier light with the test-mass mir- 
rors. A differential (common) arm-length change induces 
phase-modulation fields that only emerge at the dark 
(bright) port. Correspondingly, fluctuating modulation 
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fields that enter the interferometer from the dark (bright) 
port only interact with test-mass mirror displacements in 
the differential (common) mode. Homodyne detections 
are made at both, the dark and the bright port, each 
with a certain frequency-independent quadrature phase. 
In this way, we have two independent measurement pro- 
cesses in our interferometer, one for = {x° — i") and 
the other for x'^ = {x° +x'^). Defining p'^''^ = ztp"-)/2, 
we have = ih. As it is well known, homo- 

dyne detections of out-going modulation fields collapse 
the quantum state of the corresponding mode of test- 
mass mirror motion. In absence of classical noise, each 
mode will reach a stochastic Gaussian pure state, with 
first-order moments (i) and (p) undergoing random walk, 
and second-order moments Vxx, Vxp and Vpp remaining 
constant and minimum-Heisenberg-limited, as required 
by purity. At any instant, values of {x{t)) and {p{t)) are 
determined by measurement results in the past. In this 
sense, this Gaussian state is called a posterior state [ist . 
As long as the measurement processes for the common 
and differential modes are different in terms of e.g. sig- 
nal storage time and/or homodyne phase, the wave func- 
tions '4>'^ and il)'^ will be different, and the joint wave func- 
tion ^'(2;^a;") ^ ^""[{x"^ + x'^) /2] V/[(a;<= - a;")/2] must be 
non-separable. However, in presence of classical noise en- 
tanglement becomes less significant or even disappears. 
Note that this is analogous to creating entanglemen t by 
overlapping two differently squeezed beams on a Bs m. 

Wiener Filtering. In the following we obtain the condi- 
tional quantum state using Wiener Filtering [31, which 
applies to stable linear systems with Gaussian noise. 
More details will be given in Ref. 17[. Let y{t) be the 
Heisenberg operator of the out-going field quadrature be- 
ing measured, and x{t) be the Heisenberg operator of 
the test-mass position. If wc want to estimate x using 
the output y, the optimal filter Kx{t) is determined by 
the condition that the quantity 



Rx{t)=x{t) 



dt'Kx{t-t') y{t') 
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must be uncorrelated with the past output, i.e. 
{Rx{t)y{t') + y{t')Rx{t)) = for aU t' < t. This leads 
to a Wiencr-Hopf equation, whose solution can be given 
in terms of the single-sided (cross-) spectral densities of 
X and y: Sxy{^), Sxi^) and Sy{ft) = Sy{ft)s*{ft) with 
Sy{fl) and its inverse being analytic in the upper- half 
complex plane. Then Kx{i^) = [Sxy/sl]_^ /sy . Note 
that [...]+ stands for taking the component of a func- 
tion whose inverse Fourier transform has support only in 
positive times. The first- and second-order moments of 
the posterior state are then given by 



dt'Kx{t^t')y{t'), (3) 

§{sx-[Sxy/s;]^[Sxy/s;]l)A^) 



respectively, where y{t) is the measurement data of y{t). 
The same approach applies to the momentum operator 
p (replace every a; by p in Eqs. (I21)^(|ll)), and we further 
find 

Vx, = l"°§ "^{S^P [S.y/^;]+ [Spy/s;]l} . (5) 

Such an essentially classical filtering is justified because 
[y{t),y{t')] — (satisfied by all out-going field quadra- 
tures) and [x{t),y{t')] ^ [p{t),y{t')] = for t > t' (due 
to causality, i.e. the out-going field does not influence 
any future test- mass observable). Furthermore, it can be 
shown that the conditional moments do not change when 
different linear feedback systems are applied - as long as 
at time t the control force only depends on {y{t') : t' < t}. 

Posterior state of a single mode and squeezing. Both, 
the common and the differential mode quantum measure- 
ment process, can each be described by a set of two 
frequency-domain Heisenberg equations. 



y(fl) = sin0ai-fcos0 0,2 + a/h{x{il) + £^x) , (6) 
x{n) = -{a di + ip)/ [m{^^ + ilrr^^ - ^la)] , (7) 

with LOm the pendulum angular frequency and 7™ the pen- 
dulum damping rate. Eq. ([6]) shows that y depends on the 
in-going amplitude/phase quadratures ai 2, the motion of 
the mirror center of mass, and sensing noise ^x- The lat- 
ter may be due to optical losses, or thermal fluctuations 
of the mirror's shape, i.e. internal thermal noise, which 
makes the mirror surface move with respect to its center 
of mass. Eq. ([7]) describes the motion of the test mass cen- 
ter of mass under radiation-pressure noise (a di), as well 
as a classical force noise $,f , acting directly on the center 
of mass of each test-mass mirror, e.g. due to seismic noise 
or suspension thermal noise. The quantity a indicates 
the measurement strength, given by ad = \/4huJoP/c for 
the differential mode and ac = 2/t^/4EoJqP/c for the 
common mode, where P is the circulating power in the 
arms and r the PR-mirror transmissivity. The angle (j) 
denotes the quadrature angle of y. In the following we 
will only consider phase-quadrature readout, i.e. i/i = 0. 

From Eqs. (H]),!!?]) and with Sa^aj — % (referring to 
vacuum input) as well as S^^ and 5*^^ , determined from 
the interferometer's noise budget, we obtain the relevant 
spectral densities: Sx = [a^ + S^p]/[m'^{{fl'^ — uj'^nY + 
7^172)]^ ^ {^a/K)Sx and Sy ^ I + {a/hf[Sx + %] 
as well as Sp = m?^'^Sx, Spy = —\mQ.Sxy and Sxp = 
imU,Sx- For simplicity, we assume S^p and S^x to have a 



white spectrum, 5*^^^ = 2TimVt\ and = 2h/(rnfl^), 
even though our approach enables us to treat non-white 
spectra. The sensitivity of the interferometer is given 
by the SQL-normalized, position-referred noise spectral 
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FIG. 2: Noise spectral densities in arbitrary units with 
JIf = 27r230Hz and fia: = 27rl270Hz, and optimized ~ 
2tv 1600 Hz, Q!^ = 2-k 184 Hz. We obtain a test-mass entangle- 
ment of Eu — 0.35. 



momentum with respect to the ground state of the pen- 
dulum (recall that fla 3> Wm)- Increasing (decreasing) 
ria, e.g. by varying the circulating power or the PR- 
mirror transmissivity, results in more (less) squeezing in 
position, and more (less) anti-squeezing in momentum - 
causing the squeezing ellipse to rotate, as shown in Fig. [31 
The posterior state can also be modified when a different 
output-quadrature phase </> is used. 
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density in the detection band, i.e. 3> uj„n^n 



S'total _ 1 

^ 2 



5'sQL 



5'quant/5'sQL 



Sforco/SSQL Ssons/SSQL 



(8) 



The quantum noise Squant is limited from below by S'sql 
with equality reached at Via = a/VmK which sets a time 
scale for the measurement process. The force and sensing 

noise 5'force and S^ens 

intersect S'sql at Op and Qx, re- 
spectively. For Qx/^F > 2, we have a non-zero frequency 
band (in between flp and Qx) in which the classical noise 
is completely below the SQL (cf. Fig. [21). Solving the 
Wiener Filtering problem, wc obtain in the limit of 7™, 

< n 

a 



Vxx = n/{V2mna)[{i + 2cl){i 



■2Cff\ (9) 



Vr 



pp 



nmrj„/\/2[(l + 2C|)3(l 
n/2[(l + 2C|)(l + 2CJ)]^/^ 



1/4 



(10) 

(11) 



for the posterior state, where we have defined C,f = 
^pj^a and = ^al^x- At the bright port the in-going 
modulation fields at frequencies 51 are usually not in vac- 
uum states, i.e. Sa^a^^Saiai > 1, due to technical laser 
noise, which adds for the common mode {Sa^ai — l)/2 
and (Sajaa — l)/2 to (p and C^j respectively. Note that 
Eqs. (|9])- pTl) are consistent with results obtained from 
SME [18|, although with a more complicated noise bud- 
get, solving the Wiener Filtering problem will be substan- 
tially easier. 

A measure of the quantum nature of a state is the quan- 
tity U = VxxVpp - V^p, with U > h^/i, the Heisenberg 
Uncertainty Principle. According to Eqs. (|9ll- ([TTl) . 

u/{h^/A) = (1 + 2C|.)(i + 20 > (1 + 2nF/nxf, (12) 



with equality achieved when fla = \/Qx^f- We see that 
as long as (^p and Cx arc not too big. the posterior state is 
highly squeezed in position, and highly anti-squeezed in 



FIG. 3: Squeezing ellipses of position and momentum oper- 
ators with respect to common and differential mode for the 
same parameters as in Fig.[2] All second-order moments are 
normalized to the ground state of an harmonic oscillator at 
frequency {ni+nt)/2. 

Entanglement. After obtaining the individual common- 
and differential-mode posterior states, we assemble the 
posterior state of the entire system. The combined 
(4 X 4) covariance among (a;°,p°, x",p") reads V = 

(C(4e,Kn),(Ke,K„)) with 



ce 1 



(K^ + Vxi)/^ {Vx^p + V,%)/2 \ 

\ xp ' xp } I pp ' pp / 

{yxx~yxx)i^ {yxv~K)i'^ \ ^^A^ 

\ xp xp) I pp pp / 



A computable measure of entanglement, the logarithmic 
negativity, for an arbitrary bipartite system was intro- 
duced in [101 . For our state 



Ej^ = max[0, — logj 2cr 



(15) 



where cr" = y (S - - 4dcty)/2 and E = dot Kn + 
det yoo-2det Ko- If we insert Eqs. (O-ini) into Eq. p5l) . 
only three parameter ratios remain: Vl'^^'^ /Vlp for common 
and differential mode, respectively, and ^xI^f- The lat- 
ter turns out to be the crucial factor for the existence of 
test-mass entanglement. The parameters in Fig. [21 result 
in £V — 0.35. Recall that there exists a frequency win- 
dow with sub-SQL classical noise iff ^x/^f > 2. How- 
ever, the existence of entanglement sets a slightly higher 
threshold value for this frequency ratio. In Fig.[4l we 
plot the maximum achievable Ej^ (varying fl"^ and f2J^) 
as a function of ^x/^f for different strengths of techni- 
cal laser noise. Remarkably, even for a value of 10 dB 
above vacuum noise, the generation of entanglement is 
still possible. 

Feasibility of achieving sub-SQL classical noise. Cur- 
rent experiments 0, B H [lot aim to reach or even 
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FIG. 4: Logarithmic negativity versus Qx/i^p, maximized 
with respect to Q% and 0,^ (example values marked) for no 
technical laser noise (solid line), 5dB (dashed line) and 10 dB 
(dotted line) technical laser noise, in both the amplitude and 
phase quadratures. 



surpass the SQL by reducing the classical noise around 
rifj which is in fact not the optimum strategy for strong 
test- mass entanglement (cf. Fig. [5]). Furthermore, these 
experiments only consider a single readout. Nevertheless, 
they provide valuable information whether sub-SQL clas- 
sical noise is achievable at all. Regarding the force noise, 
it is feasible to suppress the dominant suspension ther- 
mal noise below the SQL in the detection band even at 
room temperature. This can be realized by a low pen- 
dulum frequency with a high quality factor tOm/^m- In 
Advanced LIGO [1], we expect ^F/{2Tr) ~ 40 Hz and 
ria/ (2Tr) marginally higher. In lab-scale experiments 
il.F/{2iT) = 400 Hz and fla/{2TT) of several kHz have al- 
ready been realized Regarding the sensing noise, 
optical losses need to be discussed In our setup 
the main contribution arises from photodetection loss be- 
cause no high-finesse cavities are involved. Assuming a 
total optical loss in power of e = 2%, this on its own 
would imply il^ = 14f2ct, which is more than sufRcicnt. 
More challenging is mirror internal thermal noise. In 
Advanced LIGO, we expect ^Ix to roughly coincide with 
ila- However, cooled test masses will allow to suppress 
this noise and increase Qx- In the CLIO interferometer 
currently being commissioned, the theoretical thermal- 
noise budget suggests flx/^F ^ 10 at a temperature of 
20 K . In the configuration outlined in Ref . [l^ , we 
have f2a,/(27r)=50Hz, which is a factor of 8 below ^Ip. 
With a new test-mass with lower mechanical coating loss, 
we expect fix to increase by a factor of ~ 1.5. Increas- 
ing the beam spot size by factor of 3, will provide an 
additional factor of 3^ [lO|, allowing Qx/flp ^ 1-7, even 
at room temperature. Other strategies are currently be- 
ing investigated to further suppress the internal thermal 
noise Laser amplitude noise can be suppressed suf- 
ficiently with current technology 



ity has already been reached for the injected carrier light 
of current GW detectors. 

Conclusion. In this Letter, we obtain a general and 
quantitative condition for the possibility of generating 
macroscopic entanglement between two test-mass mirrors 
of a Michelson interferometer - in terms of the interferom- 
eter's classical noise. We have found that these classical 
noise sources are required to surpass the SQL only moder- 
ately, well within the reach for experiments to be carried 
out in the GW community within the next decade. 
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Regarding laser 
phase noise, a level of 10 dB above vacuum noise - for 
an input laser power of e.g. 1 W at a measurement fre- 
quency of 1 kHz - corresponds to a frequency noise of 
about 10~^Hz/VHz. Ref. [ill suggests that this stabil- 
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